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Abstract 

A method for calculating coefficient functions of the operator prod- 
uct expansion, which was previously derived for the nonsinglet case, is 
generalized for the singlet coefficient functions. The resulting formula 
defines coefficient functions entirely in terms of corresponding sin- 
glet composite operators without applying to elementary (quark and 
gluon) fields. Both "diagonal" and "nondiagonal" gluon coefficient 
functions in the product expansion of two electromagnetic currents 
are calculated in QCD. Their renormalization properties are studied. 
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1 Introduction 



The light-cone (LC) operator product expansion (OPE) [T] (see also [2j [3]) 
permanently receives much attention, since it enables us to separate contri- 
butions to cross sections coming from large and small distances in a variety 
of processes. It was proposed as a generalization of the OPE at short dis- 
tances [1] in order to describe deep inelastic scattering (DIS) of leptons off 
nucleons. 

In Refs. [5] the T-product of two scalar currents near the LC was defined 
in term of so-called bi-local light-ray composite fields. Later on it was shown 
that the local LC expansion can be obtained by performing a Taylor expan- 
sion of the nonlocal one [5]. The nonlocal expansion is more general, but in 
the present paper we restrict ourselves to considering standard local OPE. 

The OPE for the T-product of two electromagnetic currents is of partic- 
ular importance to practical application. It can be written in the form (see, 
for instance, [7]): 
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where the dots denote contributions from other Lorentz structures. This ex- 
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composite operators and singlet (quark 0™ 

posite operators. The quantities C^f, C^ ll and C^ nl are called coefficient 
functions (CFs) of the OPE. 

As a rule, the internal sums in I are omitted in (JT|) since neglected terms 
do not contribute to DIS structure functions (see, for instance, [3]). Below 
we will discuss this point in more details. 

The standard approach to calculations of the OPE CFs is to apply for 
perturbation theory by considering the scattering of leptons off elementary 
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(quark and gluon) off-shell fields. In Ref. [7| a new method for calculating 
CFs was proposed which does not explicitly depend on elementary fields, but 
instead defines the CFs entirely in terms of Green functions of the currents 
and/or composite operators. 

In our previous paper j7] the nonsinglet case was studied. In the present 
paper we generalize our results for the singlet case. In Section [2] we derive 
a closed representation for the singlet CFs in terms of vacuum matrix ele- 
ments of the composite operators. In Section [3] we calculate the singlet CFs 
in perturbative QCD and demonstrate that our main formulae not only re- 
produces well-known expression for the gluon CF, but enables us to obtain 
the CFs of all gradient singlet operators in the OPE. The renormalization 
of singlet quark and gluon composite operators and their CFs is considered 
in Section HI A number of useful mathematical formulae is collected in Ap- 
pendix A (integrals) and Appendix B (sums). 

2 OPE coefficient functions and vacuum ma- 
trix elements of composite operators 

The cross section of deep inelastic lepton-nucleon scattering (DIS) is related 
with the hadronic tensor (see, for instance, [3]) 

W^ip, q) = 2n 2 Jd 4 xe i « x (p\TJ™(x)Jr(0)\p) • (2) 

Here \p) means a nucleon state, and J^ m (x) is an electromagnetic current: 

J7(x) = W(x)7 M <S*(s) , (3) 
where ty(x) is a quark field. The electric charge operator in (J3|), 

Q = i(A 3 + -±=A 8 ), (4) 

obeys the equation 

were A a (a = 1, 2, . . . 8) are the Gell-Mann matrices, Sp(A a ) = 0, Sp(A a A fe ) = 
25 a b, and A is the identity matrix. 
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For DIS of a charged lepton, the hadronic tensor (J2J) has two independent 
tensor structures [3]: 

W, v (p, q) = (-g, v + q -^jf\ -L Fx(xb, Q 2 ) 

+ (ph ~ \ P» - Qv^f) f ^ x bi Q 2 ) > ( 6 ) 

where Q 2 = —q 2 , 

x B = Q 2 /2pq (7) 

is the Bjorken variable, and structure functions Fi, F<i depend on these in- 
variant variables. In the Bjorken limit, Q 2 — > oo, xb is fixed, the structure 
functions Fi^{Q 2 , x) are defined via one-nucleon matrix elements of the com- 
posite operators which enter OPE (fl|). 

Near the light-cone, leading contributions to matrix elements come from 
twist- 2 operators. In QCD, the nonsinglet quark twist-2 (traceless) gauge- 
invariant operators^ are of the form (1 ^ I ^ m): 

O a Ns,'L^ = im ~ 1S 9 ^ ■ ■ ■ 9^)1^ . . . D^(x) 

+ (terms proportional to g^^) , (8) 

where operator S means a complete symmetrization in Lorentz indices, 

D lt = d li + igt a Al (9) 

is a covariant derivative, and A^(x) is a gluon field. The singlet quark twist-2 
operators (1 ^ I ^ m), 

OtL:M = ^~ 1S d ^ ■ ■ ■ d^ixh^D^ . . . D^{x) 

+ (terms proportional to <7 MiM ) , (10) 

can mix with the gluon twist-2 operators (1 ^ I ^ m — 1) 

OvL..,M = * m ~ 2s s p^ +1 • • • • • • 

+ (terms proportional to <7 WM .) • (11) 

1 Non-gauge-invariant composite operators in the OPE will be discussed in the end of 
section 01 
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Feynman rules for these composite operators which will be used for our fur- 
ther calculations are presented in Figs. [TH3J They has to be considered as a 
generalization of well-known Feynman rules j8] for the case p ^ 0. 
If the OPE ([T]) is applied to DIS, only operators of the type 

0^8 tlkL ..^Jx) = i^SVix^D^ . . . D IXm X an if(x) 

+ (terms proportional to g^ iN ) , (12) 

Of^^Jx) = i^S^x^D^ . ..D^(x) 

+ (terms proportional to , (13) 

O^^Jx) = ^- 2 SSpF WQ (x)^ 2 . ..D^F-Jx) 

+ (terms proportional to g MiAlj ) (14) 

are relevant. It is due to the fact that any composite operator O™' ^ m 
(A = NS, F, V) with at least one full derivative^ gives no contribution to a 
forward matrix element (p\ O™' ^ m \p)- In our notation, 

U NS,ti 1 ...tJ, m ~ U NS,m...^ m ' 
ryn fY 71 ' m 

rym s\m,m— 1 f-IK\ 

U V, l X 1 ... l Jt m - V, fjLl...fjL m ■ l i0 J 

In what follow, the operators ( fl5l) will be called "major" or "diagonal" com- 
posite operators, while the quark operators with 1 ^ I ^ m — 1 and gluon 
operators with 1 ^ / ^ m — 2 will be reffered to as "nondiagonal" composite 
operators. Correspondingly, we define: 



'm ^771,771 > 

c F = C F 

m m,m ' 



For nonforward matrix elements (for instance, describing deeply virtual 
Compton scattering), all composite operators contribute. Namely, we have 



2 As one can see from Eqs. (0, (Ti0|) . (fTTj) . the total number of full derivatives is equal 
to m — I or m — I — 1 for the quark or gluon composite operator, respectively. 
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the relations]^ 

(p + A| O^..,,,. \ P ) = A w+1 . . . A, m (p + A| O a / s ^ \p) , 

(p + A l °fL-Jp) = A « +1 • • • v. (p + A i > 

(P + A l C^...^ [p> = A w+1 . . . A^ (p + A| 0^... w [p> . (17) 

The "major" operators in the RHS of Eq. (ITTj) are defined above in Eq. (j!5|) . 

As usual, we assume that C^ ; (x 2 ) are tempered generalized functions 
(this is explicit in perturbative calculations), so the symbolic relation 

^ . . . = (-2,)- i^^- (-g 2 )™ (_) (18) 

holds in connection with the Fourier transform in (pQ). 

The approach, developed in our previous paper j?j for the nonsinglet CFs, 
should be generalized for the singlet case. To do this, let us 

1. take T-product of both sides of the OPE (Q]) by a singlet composite 
operator 0^ k v n ( z )i with A = F or V, 

2. imbed all resulting operator products between vacuum states. 

As a result, we obtain from Eq. (CQ) the following relation between vacuum 
matrix elements of the operator products and OPE coefficient functions: 

s og m 

< m =0 1=1 \ y / 

oo m l .. 
_ 1 n^' m ~ 

m=0 1=1 \ 1 J 

x J dSe^TO^^JO) O^...^)) + • • ■ , (19) 

3 Here (p | and |p + A) are one-particle states with 4-momenta p M and (p + A) M . 
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where C^^q 2 ) is a Fourier transform of '• 



/ m,l 
m 



CiM 2 ) = ^ (-<? 2 ) m (JpY Jfa^Chtf) , (20) 
and a new notation, 

J™{x) = W(x) 7m A°*(x) , (21) 

is introduced. In other words, only a singlet part of the product of two 
electromagnetic currents (see Eqs. (jSJ)) gives a contribution to (|T9|) . 
Let n M be a light-cone 4- vector not orthogonal to 4- momentum p^: 

nl = 0, pn^O. (22) 

Throughout the paper, we will work in the limit 

p -> , p 2 < . (23) 

Let us underline that the limit p 2 — > does not assume the limit p^ — > 0. On 
the contrary, given = %—p^ q 2 /[pq + \J (pq) 2 — p 2 q 2 ], one gets pn ~ pq at 
p 2 ~ 0. 

It is useful to convolute vacuum matrix elements with the projector 

(24) 



(pn) n ' 

and define the following invariant structure, 



n 



" {pn) : J d^xe^ J d^e^(Tj; m (x)Jr(0)Of Ui ,, Mn (z^ 

= ~g»vF2> k (u,Q 2 ,p 2 ) + .... (25) 

It depends on invariant variables p 2 , Q 2 and dimensionless variable 

u = l/x B = 2pq/Q 2 . (26) 

The vacuum matrix element of the T-product of two composite operators 
has the following Lorentz structure (A, B = F, V): 



Jd l ze*(TO$^{0)O'fc l ^(z)) 



= 2 Pm ■ ■ ■Pv.mPvx ■ ■ -PvAOb' O n A )(P 2 ) 

+ (terms proportional to g^p 2 , g„ lUj p 2 , g^p 2 ) ■ (27) 
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Equation (|2"7|) means: 



(pn) 7 



Note that both F2' k (cj , Q 2 , p 2 ) and (0£''O^' fe )(p 2 ) are dimensionless quanti- 
ties. 

Let us note that at p 2 — > vacuum matrix elements of composite op- 
erators of higher twists are suppressed by powers of p 2 with respect to the 
vacuum matrix elements of twist- 2 operators (I28I) . Thus, our approach en- 
ables us to isolate a contribution from twist-2 operators. 



2 P » 1 ...p» m (0™> l O n /)(p 2 ). 



(28) 



At fixed Q 2 and p 2 , 3-point Green function (TJ^ m J^ m 0^ k ) has a discon- 
tinuity in the variable (q + p) 2 for (q +p) 2 ^ (that is, for to ^ 1). By using 



the dispersion relation for F2 ,k (cu,Q 2 ,p 
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V 

= — V w m / du'co 
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one can derive from Eqs. (|T9l and (1251) . (!28|) : 

E^w 2 /^)(o^ fc )( P 2 /^ 2 ) 

m— 1 

+ E^W 2 //i 2 )(o^of)( P 2 //i 2 



, ImF/ (w ,g ,p 



i, ,/— m— 1 



disc^'V^P 2 ), (29) 



z=i 



J p 2 -s-0 



1 

^- y dxsx^' 1 disomy F^ k {x B ,Q 2 /p 2 ,p 2 /^ 2 ) 



(30) 



Strictly speaking, possible divergencies must be subtracted from the dis- 
persion relation f l29|) . However, it does not alter our scheme provided the 
integrals in the r.h.s. of Eq. (130]) converge (remember that m > 2). 

In ( |30l) we took into account that both matrix elements of renormalized 
composite operators and CFs depend on the renormalization scale /i. In 
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what follows, we take \i to be equal to the regularization scale /x, which arises 
in dimensional regularization [9], when one changes an integration volume, 
d 4 k fL^- D ^d D k. 

Both sides of Eq. (130]) has no dependence on n except for the trivial 
factor (—1)™- By setting k — 1, 2, . . . 2m — 1, we thus obtain a set of 2m — 1 
algebraic equations for the singlet OPE CFs C^ l (1 < I < m) and C^ nl 
( 1 < / < m 1 ) o 

Formula ( 130]) gives an operator definition of the OPE CFs in term of 
vacuum matrix elements of composite operators^ It is important to stress 
that our definition of the OPE CFs is unambiguous and it does not lean 
on a notion of quark and gluon distributions. The latter are defined via 
nucleon matrix elements of the quark or gluon composite operator, while the 
coefficient functions are independently defined via vacuum matrix elements 
of the product of composite operators. 



3 Calculations of singlet coefficient functions 
in perturbative QCD 

The formula (130]) is a generalization of a corresponding formula for a nonsiglet 
case which was derived in our previous paper [TJ: 



p 2 -s.O 



i 

l~ J dx B x^~ l disc (p+q) 2F%s(xB,Q 2 /p 2 ,p 2 /n 2 ) 



(31) 



J p 2 ->0 



By using this formula, the following expressions for the nonsinglet CFs were 
calculated in QCD [7J: 



C 
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(32) 



4 Since n > k, the index n must be chosen larger than 2m — 1. 

5 The electromagnetic current is a particular case of a quark composite operator 



with zero anomalous dimension. 
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and 
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(33) 



for I = 0, 1, ... , m — lj§ Here and in what follows superscript "(0)" means 
that a corresponding quantity is calculated in zero order in strong coupling. 
In the next order in a s , we obtained the following expressions [TJ: 
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(35) 



for Z = 0,l,...,m — 1. Here and in what follows superscript "(1)" means that 
a corresponding quantity is calculated in the first order in strong coupling 
constant. 

Let us stress that we did not demand from the very beginning that the 
"major" CF, C^ s m , should be equal to zero for odd m (see Eqs. f[3"2"j) and (13"4"|) ). 
On the contrary, it is a consequence of the fact that electromagnetic interac- 
tions conserve P-parity. Remember that DIS structure function F 2 (xb,Q 2 ) is 
an even function of Bjorken variable xb, and its nonzero moments, F%{n, Q 2 ), 
are defined by quantities C^ n (Q 2 //J 2 )(p \ 0^ n \p)(fx 2 ), with even n. That is 



'Everywhere (") denotes a binomial coefficient. 
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why we expect that the gluon "major" CF should be also proportional to the 
factor [1 + (-1)™] (see flST} below). 

For a convenience, let us for a moment rewrite our main relation ( I3T]) in 
simbolic form: 

{J JO A ) = C F (0 F A ) + C v (O v O A ) . (36) 
Then we get from (136]) : 

[(JJOv)}^ = [C F f ] [{0 F O v )f ] + [C v } {0) (OyOv)]® . (37) 
Since [(O v O F }} (0) = [{JJO V }] (0) = 0, while [(O v O v )]^ is nonzero, we get: 



m— 1 

E|<5, 

1=1 
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J m,l 
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(0) 



0. 



(38) 



Equality f[38|) is valid for all integer m, n and 1 ^ k ^ n — 1. Thus, we 
conclude that 



l(o) 







(39) 



(40) 
(41) 



for all integer m, and 1 ^ I ^ m — 1. 
Analogously, we obtain from ( 136]) : 

[<JJO F >] (0) = [^] {0) [<O F O F >] (0) + [C^] (0) (O v O F )] i0) , 
[(jjO F )] (1) = [c F ] {0) i(0 F F )] w + (O F O F )] (0) 
+ [C v ] {0) [(O v O F }} {1) + [C V } W (O v O F }} {0) . 
Taking into account that [C v } {0) = (J3H]) and [(O y O F )] (0) = 0, we find: 

[(JJO F )] (0) = [^] (0) [(O F O F )] (0) , (42) 

[(JJO F )] (1) = [C F } {0) [(0 F F )] W + [C*] (1) (O F O F )] (0) . (43) 

These equations are identical to those derived for the nonsinglet quark CF 
in our paper [7]. As a result, we find that the singlet quark CFs coincide 
with the corresponding nonsinglet quark CFs in zero and first order in a s ^\ 



^111,1 


(0) 




(0) 




^111,1 


(1) 


r>Ns 


(1) 


(44) 



7 Note, however, that [C* ,] (n) ^ [C^f] (n) , for n ^ 2 
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with [C^f] (0) and [C^f] (1) given by expressions (ESD- 

Now let us turn to QCD calculations of the gluon CFs in the first order 
in strong coupling constant by using our main formula fl30l) . We work in 
the dimensional regularization [9J and use the MS-scheme to renormalize 
ultraviolet divergences. All results of our calculations are gauge invariant 
since we sum all diagrams in each order of perturbation theory Let us 
remember that in order to find the OPE CFs, we have to retain only leading 
terms in the limit p 2 — > 0. This significantly simplifies the calculations. We 
will restrict ourselves by considering leading terms in ln(Q 2 /fi 2 ), although 
our main formula (|30|) enables one to calculate subleading terms as well. In 
other words, along with the limit p 2 — > 0, we are interested in large values of 
variable Q 2 . 

Starting from ( 136]) . we can schematically write: 



[(JJO V )] W = [C F } (0) [(0 F V )] W + [C V f } (O V O V )] {0) 

In full detail, Eq. ( 145]) looks like the following: 



(45) 
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(p+q) 2P V 
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(46) 

p 2 ->0 

. (jSSD), while the 



The quantities {C^^ are already known (see 
other terms in (T4"6"|) should be calculated. 

The propagator of the gluon composite operator is shown in Fig. HI and 

one gets 



-1) 



n+l 



In 



—p* 



B(fc + 2,Z + 2). (47) 



The vacuum matrix element of the product of two singlet composite op- 
erators is given by the diagram in Fig. |5j The calculations result in the 



^Everywhere B(a;,y) means beta-function. 
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following expression: 

(OpoF)Mtf/iM 2 )=i(-l) n C I 



47T 3 
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k k + 1 k+2 



k + l + 2 



B(Jfe + 3,0 



k + l + 1 



-B{k + 2,l) 



2 1 

+ 



/ - 1 



k k+1 k + 2 J 1(1 + 1) 
Now we are able to calculate the second term in Eq. ( 1461) : 

(0) 



E[4[ (°W> (1) (pV; 
i=i 



i(-D n C f 



32tt 3 
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— p 2 /J m(m + l)(m + 2) 



x 



m — 2 



2{m + 2) ^m{m + l) + 2 m[m + 2) 



k (k + l)(k + 2) k + m + 2 k + m + 1 
[(m - 2)B(ife, m + 3) + m(m + 2)B(k + l,m + 2)] 



(48) 



(49) 



Summation in I was made with the help of Eqs. (IB. lj) - (IB. 81) from Appendix B. 
The relation between beta-functions (integer k, m), 



2 2 1 \ 1 

+ ] B(k + 3, m) - - B(k + 2, m) 



1 



k k + 1 ' k + 2) v " ' " ' k K " ' ' m(m + l)(m + 2) 
x [(m - 2)B(k, m + 3) + m(m + 2)B(fc + 1, m + 2)] , (50) 



was also used. 

The diagrams which contribute to vacuum matrix element with two cur- 
rents are shown in Fig. [6j Omitting details of calculations, let us give the 
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result: 



J fa" 1 ' 1 disc (p+9)2 (J J Oy k )( 1 \x B ,Q 2 /p 2 ,p 2 / H 2 ) 





In (-^l ln^ 
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—p 2 J m(m + l)(m + 2) 
m(m + l) + 2 m(m + 2) 



fc ' (A; + l)(fc + 2) ' A; + m + 2 £; + m + 1 
- [(m - 2)B(fc, m + 3) + m(m + 2)B(fc + 1, m + 2)] I , 



(51) 



In order to obtain Eqs. ( l4T|) -( l5"Tj) . we used integrals (IA.lj) - (lA.3j) from Ap- 
pendix A. 

Equations (14"5]) -(]5T ]) result in a set of equations for the gluon CFs. Namely, 
for any integer m ^ 2 we obtain algebraic equations for [C^ i(Q 2 / fJ 2 )]^, with 
1 < I < m- 1: 



m— 1 



2=1 



X 



(1) B(A; + 2,/ + 2) = C F ^ln (%) — L 

47r \p / m(m + l)(m + 2) 



m-2^ 2(m + 2) m(m + 1) + 2 m(m + 2) 



k (k + l)(k + 2) k + m + 2 k + m + l_ 



- [(m - 2)B(fc, m + 3) + m(m + 2)B(fc + 1, m + 2)] 



(52) 



Note that these equations holds for all integer k ^ 1, but for our purposes it is 
enough to consider only m—1 equations corresponding to& = l,2,...m— l|j 
The solution of equations (152]) is a sum of two terms one of which is 
nonzero only for even m, while another is nonzero only for odd m: 



clAQ 2 /^ 2 ) 



a) [i + (-i) m ] 



2 

[i-f-r 



i(i) 

even 
(1) 

odd 



(53) 



9 The other equations which correspond to k m, will be also satisfied, as one could 
see from an explicit expression for our solution (|53p - (|55|) . 
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The first term in (1531) is defined for 1 ^ Z ^ m — 1 



(1) 



x 



1 



2 2 

+ 



m m+ 1 m + 2 



47r y /i 2 
-1)' 



m 
Z + l 



— (m — I) 



m + 1 



:-D'['™ »') -„»-/-i. 



(54) 



while the second term in ([53]) is nonzero for 1 < Z < m — 2: 



l(i) 



a 



C F — In 



odd 



47T 



/x 2 / m(m+l)(m + 2) 



x J (-1)' [(/ - m + l)(m + 2) + rn(m + 1) + 2] 



[(/ - m + \){m - 2) + m(m + 1) + 2] 



m 
Z + l 



(55) 



Note that (C^Jodd (1531) gives no contribution to C^ ll ( {53"]) for Z = m — 1 
due to relation [(— l) m — 1][(— l) m + 1] = 0. It is rather easy to demonstrate 
that (|53"1) does obey set of equations (152]) for any m ^ 2, k ^ 1, if one uses 
formulae (1B.9I) - (1B.12I) from Appendix B. Indeed, these formulae lead us to 
the relations: 

m— 1 m-1 

£ B(fc + 2, Z + 2) = £ fc( W)l „„ B ( fc + 2, Z + 2) 

L J even L J odd 

Z=l «=1 



C F — In ^~ 



47T 



1 



/x 2 / m(m + l)(m + 2) 



m - 2 2(m + 2) m(m + 1) + 2 m(m + 2) 
fc (fc + l)(fc + 2) fc + m + 2 £; + m + 1 



[(m - 2)B(fc, m + 3) + m(m + 2)B(fc + 1, m + 2)] } , 



(56) 
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In particular, it follows from Eq. ([53]) . (J5IJ) that 



(1) 



47T 







1 

m 



+ 



771+1 m + 2 



(57) 



As one can see, "major" CF (jo7|) is defined by well-known anomalous dimen- 
sion [10@ 

or, „ (m + 1) 2 + (m + 1) + 2 



7™v 



■cv 



(58) 



7r m(m + l)(m + 2) ' 

Thus, we have reproduced the standard expression for the "major" CF, 
m -i(Q 2 f^ 2 ), and, which is more important, have calculated "gradient" 



gluon CFs, C^CQVa* 2 ) (/ 
coupling a s . 



1,2. . . . m — 2), in the first order in strong 



4 Renormalization of singlet composite oper- 
ators and coefficient functions 

Let us consider in more detail products of renormalized composite operators 
and corresponding renormalized CFs which enter the OPE of two electromag- 
netic currents flTJ). Both singlet composite operators, O™' and O™' 1 , depend 
on the renormalization scale /io, and mix with each other under rescaling 

fl -> \L\ 

1=1 1=1 

k+1 k 

i=i i=i 

(59) 

10 In order to obtain the expression for jpy in standard notations, one has to replace 
index m by (n — 1) in (|58p. 
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where Z are the matrices of a finite renormalization of the composite opera- 
tors. In its turn, (159 p means that 



and 



k=l+l 

m 



(60) 



(61) 



k=l 



for / = 1,2,... m- 2. 

Since the quantity /x is an arbitrary scale, one can put //q = Q 2 in (160 
flBTj) . and obtain: 



and 



Cm,m-l(Q 2 1 ^) — Cm,m(l) )m-l {Of I ^) 

+ CV{1){Z VV )2Z\{Q 2 /^) 



/c=«+i 

m— 1 

+ E0 1 ) (^ yy )f(W), 



(62) 



(63) 



for / = 1, 2, . . . , m — 2. As a result, we find equations for the leading parts of 
the gluon CFs in the first order in the strong coupling (1 ^ / ^ in — 1): 



clAQ 2 /^ 2 ) 



(i) 



(0) 

E [(^ y )f(W)_ 



(i) 



(64) 



fc=Z+l 



In deriving relation ( 1641) . we took into account that K ( ]^= for all m and 
1 ^ I (1591). By using explicit form of [C^ )fc ] (0) (see flU, flH2])), these 
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equations can be written as follows (1 ^ I ^ m — 1): 



(1) 



[i+(-i)H p^)r(Q7^) 

m— 1 



(1) 



fc=/+l 



(i) 



(65) 



Note that the last term in fl65l) is identically zero at I — m — 1. 

The mixing of singlet quark operators ffTUl) an d g luon operators ffTTj) is 
defined by the set of diagrams presented in Fig. [Tq The sum of divergent 
parts of these diagrams is given by the expression: 

\~9iiv kn(k + p)n - n^n u k(k + p) + + + A^i, fen] 

2 \e fc-i 



x C, 



X 



a s l f fi 



57T £ 
1 

k ~ k 



-p- 



+ 



1=0 

2 



n-l-l 



1 ' A; + 2 
1 2 



(kn) (pn) 
1 



n-J-1 



k+1 k+2 



k + 1 

k 

l + l 



(k-l-1) 



k + 1 V l + l 



k-1 



(66) 



In deriving (ESJ), basic integrals (1A.4[) - (1A.7|) from Appendix A were used. Let 
us remember Feynman rule for the unrenormalized gluon operator Oy k (see 
Fig. ED: 

(-l) n - k ~\kn) k -\pn) n - k - 1 [ - g^ knik + p)n - n^n u kik + p) 

+ n^k u (k + p)n + k^n u kn] . (67) 

As a result, the matrix of the finite renormalization in Eq. (161| has the 
following form (1 ^ / ^ k — 1): 



(Z^)f(QV^) 



(i) 



+ 



fc k+1 k+2 



-(-1) 



1 



(*-J) 
2 



t< 2 



k k+1 k+2 



k+1 

k 

l + l 



(k-l-1) 



k 



k + l\l + l 



(68) 



xl The expressions for unrenormalized singlet quark and gluon composite operators are 
presented in Fig. [1] and Fig. [2j respectively. 
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In particular, it follows from (j68p : 



>FV\k 



(i) 



+ 



k + 1 fc + 2 



(69) 



As one can see from (j5"9"|) . [(-^ Fy )fc_i] < ' 1 ' ) = 0, for = 1. It is a consequence 
of the fact that the singlet operator \l/(a;)7 At \l/(x) do not mix with the gluon 
operators!^ 

From ( 16"5]) . ( 161?]) the expression for the "major" gluon CF (1571) follows 
which was obtained in the previous section by solving set of equations (152]) . 
In order to obtain [C% t (Q 2 / fJ?)]™ for 1 ^ / ^ m — 2, one should calculate 
the sum in k in Eq. (165]) . It can be done with the help of formulae ( ]B.13j) - 
( 1B.19P from Appendix B. As a result, we come to expressions (153"]) . (151j) . ( 155]) 
derived above on the basis of our main formula (150]) . 

Finally, by using Eqs. (IB.9[) -f lB. 12|) from Apendix B, one can obtain: 



TO— 1 



(0) 



■^//•\'V"iV / , i 



(1) 



(Q 2 //2 2 ) B(k + 2,l + 2) 



47T 



m(m + l)(m + 2) 



+ 



2(m + 2) 



fc ' (Jfe + l)(Jfe + 2) 
m(m + 1) + 2 m(m + 2) 



+ m + 2 + m + 1 

[(m - 2)B(fc, m + 3) + m(m + 2)B(fc + 1, m + 2)] 



(70) 



Correspondingly, with the help of formulae (lB~T3]) - (lB~T9]) and dB~9l- (lB~T2]) . 



12 



Obviously, it should takes place in all orders in a a 
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we are able to find: 



m— 1 m—1 

S [^™> fc 



1(0) 



=1 fc=Z+l 



z^)f(g 2 /^)] (1) B(fc + 2,/ + 2) 
Q 2 ~ 



777(777 + l)(m + 2) 



m — 2 



2(m + 2) 



fc (Jfe + l)(Jfe + 2) 
777(777 + 1) + 2 777(777 + 2) 

k + 777. + 2 + 777 + 1 

[(m - 2)B(fc, 777 + 3) + 777(777 + 2)B(fc + 1, m + 2)] 



+ 



(71) 



Thus, we have successfully reproduced equations ( 154)) and (|55|) . 

It is well known that any Green function with an insertion of one com- 
posite operator is multiplicatively renormalized [TT], while Green functions 
with insertion of two ( or more ) composite operators need additive countert- 
erms [12]. Nevertheless, as was shown in [7], renormalization group equations 
for the CFs have no additive terms, provided the corresponding composite 
operators have zero vevf^l 

It was found that some gauge-invariant singlet composite operators can 
mix with gauge-variant ones under renormalization [T3]. This problem is 
present for the simplest of these operators, the energy-momentum tensor 
M „, already in the leading order in strong coupling. 

Let O and iV represent a set of gauge-independent and non-gauge- indepen- 
dent operators, respectively. It was proven that renormalized and unrenor- 
malized operators of these types are related by a triangular matrix 



R 



Zoo 




(72) 



u 



where Zab are matrices. In other words, gauge-variant operators do not 
mix with gauge- invariant operators under the renormalization. Correspond- 
ingly, anomalous dimensions of gauge-independent operators can be deter- 
mined from matrix Zqo alone. Moreover, in so-called physical (axial) gauge, 



13 The renormalization properties of the composite operators with nonzero vev were 
studied in [T3l. 
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n^A^ = 0, n M being a constant light-like vector [T5], a renormalization pro- 
cedure does not require gauge-variant counterterms for the gauge- invariant 
composite operators at all [16] . 

It follows from f J72|) that in the OPE renormalized CFs of gauge- independent 
operators change under rescaling of renormalization mass, yUo — > as follows: 

Co^ 2 ) = Co(^)Zoo(^/^), (73) 

where Zoo is a matrix of finite renormalization. All physical matrix elements 
of gauge- variant operators vanish [H]@ Using a complete set of hadronic 
states |n), we find that 

(NO) = ^2(N\n)(n\O) = 0. (74) 

n 

Thus, a presence of gauge-variant composite operator in the OPE ([TJ) have 
no influence on our method of determining CFs in terms of vacuum matrix 
elements. Indeed, multiplying elements of the sum ^ m [Cq O m + C™ N m ] by 
one of the gauge-invariant operators, O n , and putting them between vacuum 
states, we exclude a contribution from the gauge-variant operators to our 
main equation fl30|) . 

Taking all said above into account, we did not include gauge-dependent 
composite operators in the OPE ([TJ). 

5 Conclusions and discussions 

As it was shown in the present paper, the singlet CFs of the OPE of two 
currents can be explicitly expressed in terms of the Green functions of the 
corresponding composite operators without explicit use of the elementary 
(quark and gluon) fields. Our main equation ([30]) is a generalization of an 
analogous formula which was previously obtained for the singlet case [TJ. It 
is necessary to stress that our formula holds in any renormalization scheme 
in contrast with other prescriptions (see, for instance, [T9]). 

As an illustration of a validity of our scheme, the gluon CFs were calcu- 
lated in QCD in the first order of the strong coupling constant. It is important 

14 The calculations made in |17) contradict this result. However, it was shown by explicit 
calculations [18] that the proof in [17] breaks down, and conclusions of paper [14] remain 
true. 
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to note that both "diagonal" CFs, C , ^ im _ 1 (<5 2 //i 2 ) (15TI) . and "nondiagonal" 
CFs, C 9 m l {Q 2 /fi 2 ) (1 < I < m - 2) ([53j)-([55D, were simultaneously obtained. 
The renormalization of these composite operators and their CFs were also 
considered. 

For further discussion, let us rewrite a set of equations for the singlet 
quark and gluon CFs in the following symbolic form: 

(JJO q ) = C q (O q O q )+C g (O g O q ), 

(JJO g ) = C q (O q O g ) + C g (O g O g ). (75) 

These equations must be considered as the set of matrix equations (for sim- 
plicity, summations in I are omitted). It is also assumed that the procedure 
described in details in Section [2] is applied to all matrix elements in (175|) J^I 
Let us define "reduced" matrix elements of the quark and gluon composite 
operators in the n-th order of perturbation theory (n > 1): 

{o7o q ) {n) = (o A o q )^[(o q o q )^}-\ 

(6Zo g )W = (0 A O g )V> [{O g O g )^Y l , (76) 

where A = q, g (see Figs. [HUH])- Analogously, we can define "reduced" ma- 
trix elements which contain both composite operators and electromagnetic 
currents (n > 0): 

(JJO q )^ = (JJO q )^ [{o.o^}- 1 , 

(Jl0 g )^ = (JJO B )M [{O.O^Y 1 (77) 

(see Figs. ITU1ITT1) . Note that each of the matrix elements (OaOb)^ and 
(JJOa)^ has a divergency related with a divergency of a corresponding 
Feynman graphs as a whole J^l However, these divergences cancel in the 
reduced matrix elements (O a Ob)^ and (JJO A )^ . 
From (1751) we get 

(JJO q )^ = Cf\O q O^ + Cf>{O a O q )<® , 

(JJO g )^ = Cf\O q 9 )^ + Cf\O g 9 )^ . (78) 



15 See derivation of formula ([30]) . 

16 Note that zero order matrix elements (O q (p) O q (—p)}(°\ (O g (p) O g (—p))(°\ and 
disC( p+9 )2 (J(q) J{— (p + q))O q (p))W are proportional to ln(l/p 2 ) at p 2 — > 0. 
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Since (JJO g )W = (O g O q )<® = (O q O g )<® = 0, while (O g O g )^ ^ 0, we 
immediately obtain in the leading (zero) order in the strong coupling a s : 

Cf = , (79) 
C® = (JJO q )W. (80) 

In the first order in a s , we derive the following expressions for the singlet 
quark and gluon CFs: 

cf = (jJO q )^ - (jTo q )^(o q ~d q )^) , (si) 

C« = (JJ~O g )^ - (JJO^i^d,)^ . (82) 

The first order contributions to the singlet CFs are presented in Fig. [TH and 
Fig. EU respectively. In the next order the singlet quark CF looks like 

cf = (J7o q )W - [(JJO,)® - (JJO q )QHO^V>]0&)<!> 

- [{Jjojv - (jjq;>W(p^>w]<o^>« 

- (JJO q )^ (0^ q ) {2) . (83) 
Correspondingly, the gluon CF has the form: 

Cf = (JJO g )^ - [(JJ0 9 )V> - (JJO q y°\^ q )^(^O g )W 

- [{JJo 9 ) m - (j7o q y°Ho q ^ g )W](d g ~o g )^ 

- (JJO q )M (C^Of 2 ) . (84) 

In the same way, singlet quark and gluon CFs can be calculated at any order. 

One more advantage of our approach is that it treats uniformly the CFs of 
light quarks (q = u, d, s) and CFs of heavy quarks (Q = c, b). For instance, 
singlet heavy quark CF, Cq, is given by the same formulae ( IHTj) . (!83|) . At 
the same time, the so-called double counting problem does not arise at all@ 

The leading parts of the first terms in Eqs. (1811 -(ISl). (JJC^) (n) (A = 
q, g) are proportional to ln n (Q 2 /p 2 ), while the leading parts of the matrix 
elements (OaOb)^ (A, B = q, g) are proportional to ln n (/x 2 /p 2 ). Thus, the 

17 One possible solution of double counting problem, which appears in deriving CFs of 
heavy quarks at the diagram level, was proposed in |20j . 
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role of all terms in the r.h.s. of Eqs. ( I8ip -( l84p . except the first one, is to 
cancel p 2 - dependence in the final expression for the singlet CFs: 



c 



(2) 



iV 



(2) 



hi 



p 2 



ln 2 f^ 



N {2) In 2 



P 

" 2 



-In 



N$ ] In 



P 2 J \P 2 



l< 2 



(85) 



(86) 



(n) 

where constants N A are known from explicit calculations. 

An analogy can be drawn between our formulae fl81l ) -( l84j) and diagram 
approach to calculating CFs@ Namely, the quantity (JJOa) should be asso- 
ciated with the amplitude shown in Fig. [T3J with A being the type of off-shell 
parton with 4-momentum p. The n-th order contribution to the imaginary 
part of this amplitude grows as ln n (Q 2 /p 2 ) at p 2 — >• (or, equivalently, at 
Q 2 oo). The detailed diagram analysis can be found in Ref. [21], where it 
was shown that the /independence drops out if a gauge-invariant set of QCD 
diagrams is taken into account in each order of perturbation theory. 

The quantity (OaOb) should be associated with the so-called cut ver- 
tex [22J depicted in Fig. [TSJ where upper (lower) lines in Fig. [02 correspond 
to partons of type A(B). Being integrated in 4-momentum k of the upper 
parton, this diagram in the n-th order has a singularity \n n ([i 2 /p 2 ) at p 2 — > 0. 

Our results can be applied to studying generalized parton distributions 
(GPDs)0 [23j [2H [251 12S] which permanently attract a great amount of inter- 

20| 



esto They parameterize nonperturbative parton correlation functions in the 
nucleons and interpolate between the ordinary parton distribution functions 
(PDFs), which can be measured in DIS, and the elastic form factors. GPDs 
appear in cross sections of deeply virtual Compton scattering (DVCS), hard 
leptoproduction of vector mesons, as well as in diffractive Z°-production in 
ep-collision. They were also introduced in the context of the spin structure 
of the nucleons [25] . 

If the OPE ([!]) is applied to DIS, only "diagonal" operators of the type 



^F,/^i.../i„ 



^F, fii...fi„ 



(O 



111 



O 



m,m—Y 
V, /£]... ./*n 



are important, since forward 



18 At the moment, we restrict ourselves by the "diagonal" CFs which survive in DIS. 
19 Also called off-forward or nonforward PDs. 

20 For the first time, nonforward QCD planar ladder diagrams were studied in [27] . 
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matrix elements of these operators with 1 < I < m — 1 (1 < / < m — 2) are 
zero. However, for DVCS and other processes mentioned above, all operators 
contribute proportionally to (p— ■ ■ ■ (p~P')nm- The invariant structures 
of matrix elements of these operators related to the GPDs. Thus, our scheme 
of calculating CFs of the "nondiagonal" composite operators 0™'L ^ (A = 
F, V) become quite important. 

For the first time, the very notion of nonforward distribution function 
was introduced in Ref. [28], in which the statement was made that its Fourier 
trasformation "can be interpreted as the distribution of partons in momen- 
tum fraction x and in impact parameter 6j_" . Later on, it was shown that 
GPDs in the limit when the momentum transfer is purely transverse de- 
scribe the distributions of unpolarized (polarized) partons in the transverse 
plane [29]. Impact parameter dependent PDFs satisfy positivity constraints 
which justify their physical interpretation as probability densities [29] . 

In conclusion, let us stress again that Eq. fl30|) . which defines the singlet 
CFs in term of the Green functions of the composite operators, does not 
apply to perturbation theory at all. Therefore, our results can be used for 
calculating CFs of the OPE by nonperturbative methods. 

All Feynman graphs presented in the present paper (Figs. [T ¥T5|) were 
prepared with the use of Axodraw package [30] and JaxoDraw graphical user 
interface [31]. 
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Appendix A 

In this Appendix we have presented several basic D-dimensional integrals 
(D = 4 — 2 e) which are needed for our calculations. Light-cone 4- vector 
is defined in the main text (l22p . Only divergent parts of the integrals are 
shown. 



(/' 



2\e 



d D k {kn) m i ] 

(2tt) d k 2 (k + p) 2 ~ 16?r 2 £ \-p 2 J m + 1 



fi 2 V (-1} 



{pn) 1 



(A.l) 
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(^ 2 ) 



, , dPk (kn) 1 



i 1 



{2ti) d {k + p) 2 {k-l) 2 16tt 2 s 



(_l)m 

m + 1 



p=0 



(A.2) 



(^) 



, d D k (kn) m [(k + p)n] n i 1 f fi 2 ^ £ 



(2n) D k 2 (k + p) 2 lQn 2 e\-p 2 

x B(m + 1, n + 1) , 



-l) m (pn) 



(A.3) 



where B(a;, ?/) is the beta-function, and m > 0. In the next four integrals 
m > 1 is assumed: 



(/i 2 ) e ' k^knj 



(2tt) d k 2 {k + p) 2 16ii 2 e\-p 2 J m + 2 



i 1 f fi 2 X (-I)™" 1 



m— 1 



X 



m 



—n 



M 2(m + 1) 



p 2 + p^ (pn) 



(A.4) 



I d D k k^{kn) 



i 1 



(2ir) D k 2 (k+p) 2 (k-l) 2 32ir 2 e 



-(i+ P y 



£ ( 1 Nm— 1 



1)' 



m + 1 



m— 1 



x ^(-l) p (/n) p (pn)^- 1 , (A.5) 

p=0 
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d D k k^kn) 1 



i 1 



(k + p) 2 (k-l) 2 1Qti 2 e 

m—l 



<i+pY' 



(-i) m 

m + 1 



p=0 



m—l 



l — ^T(-l) p (p + l)(p + 2)(H P (pn)^- 1 



m + 2 



p=0 



m 

III \~ — 



m—p 



p=0 



^^r(-i)w(H m ^J>, 

p=0 



(A.6) 



i 1 



cFA; k^k u (kn) m 
(2tt) d k 2 {k + p) 2 {k-l) 2 ~ 32^2 e 

f m 

x ^£(-l)W(pnr-* 



■(Z+p) J 



(m + l)(m + 2) 



p=0 



m—l 



- (Z„n„ + V*) ]T(-l) p (p + l)(Zn) p (pn)^" 1 

m—l 

+ (p^ + n^) ]T(-l) p (m - p)(Zn) p (pn)^" 1 

p=0 

m-2 

/c 2 ^(-l) p (p + l)(/n) p (pn) m ~ p - 2 

p=0 
m-2 

+ p 2 ^(-l) p (m-p- l)(Zn) p (pn) m - p ^ 2 



1 

2 



- - n M n„ 



p=0 



m-2 



+ (/ + p) 2 ^(-l) p (p + l)(m-p- 1) 

p=0 



x (Zn) p (pn) 



p (™\ m -p- 2 



(A.7) 
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Appendix B 



In this Appendix we collected formulae which are needed for calculating sums 
presented in the text and getting compact expressions Everywhere below 
(^) denotes a binomial coefficient, and ~B{x,y) beta-function. For integer 
m, n ^ 0, one has B(m + 1, n + 1) = \{m + n + 1) f m+n J] -1 . 
Let us first consider summation in index I: 

m— l / _ i\ i 

5> 1} ' (T- 1 ) B( * + 3 ' = " (_ir B( * + 3 ' m) " k + m + 2 ' (B ' 1} 

liit- 1 )' (7- /) r^r B (* + ' + 2 ) 

i=i v 7 



_(_l)"*^_B(A; + l,m + 2) 

771+1 

1 

(k + m + !)(£; + m + 2)(fc + 1) ' 



B^7:/)(7TT)W 3(M + 3) 

1 

B(A;,m + 3) 



(k + m + l)(k + m + 2)fc(fc + 1) 



-B(/c,m + 3) 



(fc + m + 2)A;(A; + l)(fc + 2; 



(B-2) 



m+l)(m + 2) 

' (B.3) 



m(m + l)(m + 3) 

' (B.4) 



21 All the formulae collected in this section were derived by using several table sums with 
binomial coefficients 132 . 
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gVl)' (7:/) "B(* + l,m) -(-I)' 



l + k 

In particular, we get from flB.5j) : 



m + k 



(B.5) 



(B.6) 



(B.7) 



.(-!)»- 



1 + 2 m + 2 K '' m(m + l)(m + 2) ' 
Next four sums in I contain the same beta-function B(k + 2, 1 + 2): 

m(m — 1) 



(B.8) 



m-l ✓ v 

E(- 1 )\; + i) B ^ + 2 ' Z + 2 ) 



(A; + m + 2)(A; + 2)(A; + 3) 



(B.9) 



£(-l)^™J(m-l-l)B(* + 2,J + 2) 



m(m — l)(m — 2) 
(k + m + l)(k + 2){k + 2>) 



m—l 



B(k + 2,l + 2) = -B(k+l,m + 2) + 



i=i 



m—l 



(k + l)(k + 2)(k + 3) 



(B.10) 



, (B.ll) 



^(/ + 2)B(fc + 2,/ + 2) = -B(jfe,m + 3) - (m + 2)B(jfe + l,m + 2) 

6 



i=i 



+ (m + l)B(fc + 2,m + 1) + 



(A; + l)(A; + 2)(A; + 3) ' 



(B.12) 
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The following four sums in k are also used during calculations. The 
Kronecker symbols 5 n>m garantee that all the sums are equal to zero for 
I — m — 1: 



m—l 



k=l+l 



m 



1 / k 



k-ljk\l + l 



■ 1 ) m k(l + l) [1 - 5o > m - l - l] > (R13) 



= -(-i) 

+ (-!)' 



m 
Z + l 
1 



1 



m + 1 m 



m(m + 1) 



m— l—l 



(B.14) 



m— 1 



= -(-i) 

+ 2(-l)' 



Z + l 
/ + 2 



1 1 , 

— <JQ,m-l-l 



m + 2 m 



m(m + l)(m + 2) ' 



(B.15) 



i-iy 



1 /m — l 



m + 1 V Z + 1 



+ (-1)- 

m 



Z + 2 



'0,m-J-l 



(B.16) 



m + 1 

Note that the latter sum is equal to zero both for I — m — 1 and / = m — 2. 



m \ 1 



m— 1 

fc=i+i 



1 

m 



:-iy( m l v ) + (-iy 



(B.17) 



m-l / \ 1 



m 



ra + 1 W + 1 



m 



m — l 
/ 



m + 1 ' 



(B.18) 
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k=l+l 



( 



k-l)k + 2 



m \ 1 



m + 2 \/ + 2 



1 fm + 1 



) 



m + 1 V + 2 



) 



+ — 




OH- 



m + 2 



1 



(B.19) 



References 

[1] Y. Frishman, Phys. Rev. Lett. 25 (1970) 966; Ann. Phys. 66 (1971) 
373; R.A. Brandt and G. Preparata, Nucl. Phys. B 27 (1971) 541; H. 
Fritzsch and M. Gell-Mann, in Broken Scale Invariance and the Light 
Cone, London, Gordon and Breach, 1971. 

[2] B. Geyer and D. Robaschik, Fort. Phys. 27 (1979) 75; A.J. Buras, Rev. 
Mod. Phys. 50 (1982) 199. 

[3] F.J. Yndurain, Quantum chromodynamics : An Introduction to the The- 
ory of Quarks and Gluons, Springer- Verlag, 1983. 

[4] K. Wilson, Phys. Rev. 179 (1969) 1499; W. Zimmermann, in Lectures 
in Elementary Particles and Field Theory, M.I.T. Press, 1970. 

[5] S.A. Anikin and O.I. Zavialov, Ann. Phys. 116 (1978) 135; O.I. Zavialov, 
Renormalized quantum field theory, Kluwer, Dordrecht, 1990; Theor. 
Math. Phys. 138 (2004) 370. 

[6] B. Geyer et al, Z. Phys. C 26 (1985) 591; J. Blumlein, B. Geyer and D. 
Robaschik, Nucl. Phys. B 560 (1999) 283. 

[7] A.V. Kisselev and V.A. Petrov, Phys. Rev. D 71 (2005) 085020. 

[8] E.G. Floratos, D.A. Ross and C.T. Sachrajda, Nucl. Phys. B 129 (1977) 
66; ibid Erratum B 139 (1978) 545; ibid B 152 (1979) 493. 

[9] G. t'Hooft and M. Veltman, Nucl. Phys. B 44 (1972) 189; C.G. Bollini 
and J.J. Giambiagi, Nuovo Cim. 12 B (1972) 20; G. Cicuta and Mon- 
taldi, Nuovo Cim. Lett. 4 (1972) 329. 

[10] D.J. Gross and F. Wilczek, Phys. Rev. D 8 (1973) 3633; H. Georgi and 
H.D. Politzer, Phys. Rev. D 9 (1974) 416. 



31 



[11] W. Zimmermann, Ann. Phys. 77 (1973) 536. 

[12] J.C. Collins, Renormalization, Cambridge Univ. Press, 1984. 

[13] CM. Shore, Nucl. Phys. B 362 (1991) 85. 

[14] S. Joglekar and B.W. Lee, Ann. Phys. 97 (1976) 160; S. Joglekar, Ann. 
Phys. 108 (1977) 233. 

[15] A. Bassetto, G. Nardelli and R. Soldati, Yang-Mills Theories in Alge- 
braic Non-Covariant Gauges: Canonical Quantization and Renormal- 
ization, World Scientific, 1991. 

[16] C. Acerbi and A. Bassetto, Phys. Rev. D 49 (1994) 1067. 

[17] R. Hamberg and W.L. van Neerven, Nucl. Phys. B 379 (1992) 143. 

[18] J.C. Collins and R.J. Scalise D 50 (1994) 4117. 

[19] K.G. Chetyrkin, S.C. Gorishny and F.V. Tkachov, Phys. Lett. 119 B 
(1982) 407; S.C. Gorishny, S.A. Larin and F.V. Tkachov Phys. Lett. 124 
B (1983) 217. 

[20] A.V. Kisselev, Phys. Rev. D 60 (1999) 74001. 

[21] D.J. Pritchard and W.J. Stirling, Nucl. Phys. B 165 (1980). 

[22] A.H. Mueller, Phys. Rev. D 18 (1978) 3705; Phys. Rep. 73 (1981) 237. 

[23] D. Muller et ai, Fortschr. Phys. 42 (1994) 101. 

[24] A.V. Radyushkin, Phys. Lett. B 380 (1996) 417; ibid B 385 (1996) 333; 
Phys. Rev. D 56 (1997) 5524. 

[25] X. Ji, Phys. Rev. Lett. 78 (1997) 610; Phys. Rev. D 55 (1997) 7114. 

[26] P. Hoodbhoy, X. Ji and W. Lu, Phys. Rev. D 59 (1998) 014013. 

[27] J. Bartels and M. Loewe, Z. Phys. C 12 (1982) 263. 

[28] A.V. Kisselev and V.A. Petrov, Singular structure of the parton distri- 
butions in semihard and soft processes, in Proc. XXIVth Recontres de 
Moriond, ed. J. Tran Than Van (Editions Frontieres, Gif-sur Yvette, 
1989), p. 213. 



32 



[29] M. Burkardt, Int. J. Mod.Phys. A 18 (2003) 173. 



[30] J.A.M. Vermaseren, Comp. Phys. Comm. 83 (1994) 45. 

[31] D. Binosi and L. Theussl, Comp. Phys. Comm. 161 (2004) 76; D. Binosi 
et al, l arXiv:0811.41131 

[32] A. P. Prudnikov, Yu.A. Brychkov and O.I. Marichev, Integrals and Se- 
ries, Vol. 1, Elementary Functions, Gordon & Breach Sci. Publ., New 
York, 1986. 



33 





Figure 1: Feynman rules for the quark composite operators 0™f Mn in the 
leading (zero) order in strong coupling a s . 
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(pn) n - k - 1 (kn) k - 1 

x [— fen (A - p)n - n^riv fe(fe - p) + n^fe,, (A; - p)n + (A; - p)^ fen] 




(-l) n - k -\pn) n - k - 1 (kn) k - 1 
x [— ff/ii/ kn(k + p)n - n M n y A; (A; + p) + n M fe„ (A; + p)n + (A; + p) M n„ fen] 

Figure 2: Feynman rules for the gluon composite operators O™^ ^ n in the 
leading (zero) order in strong coupling a s . 
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Figure 3: Feynman rules for the quark composite operators 0™f Mn in the 
first order in strong coupling a s . 
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Figure 4: The diagrams for the propagator of the gluon composite operator 
(Oy k Oy' 1 )^ in zero order in strong coupling a s . 
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Figure 5: The diagrams for the mixing of the composite operators 
(O^O™' 1 )^ in the first order in strong coupling a s . 
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Figure 8: The redefinition of the matrix element of the quark composite 
operators in the first order in strong coupling a s . 
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Figure 9: The redefinition of the matrix elements of the quark and gluon 
composite operators in the first order in strong coupling a s . 
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Figure 14: The hadronic part of the amplitude of deep inelastic lepton scat- 
tering off a parton with 4- momentum p (p 2 < 0). The dotted line means that 
the imaginary part of the amplitude should be taken. 
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Figure 15: The partonic cut vertex. The solid lines represent quarks or gluons 
fields. The "target" parton is off-shell, p 2 < 0. Integration in 4-momentum 
k is made. 
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